Generalized Berry Conjecture and mode correlations in chaotic plates 
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We consider a modification of the Berry Conjecture for eigenmode statistics in wave-bearing 
systems. The eigenmode correlator is conjectured to be proportional to the imaginary part of 
the Green's function. The generalization is applicable not only to scalar waves in the interior of 
homogeneous isotropic systems where the correlator is a Bessel function, but to arbitrary points of 
heterogeneous systems as well. In view of recent experimental measurements, expressions for the 
intensity correlator in chaotic plates are derived. 
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In 1977 Berry conjectured that the higher eigenmodes 
of a ray-chaotic Hamiltonian, in particular, a billiard, 
should be statistically indistinguishable from a superposi- 
tion of plane standing waves of all directions, with uncor- 
rected amplitudes and phases Q; the idea is also found 
in 0. One immediate consequence is that the modes are 
Gaussian random functions, with correlations given by 

(vW (x) u<"0 (x + r)) - A 2 S nm J (kr) , 

where brackets (■) represent spatial average over posi- 
tion x, A is an unimportant normalization factor, r = 
r| is a separation distance, and k is the wavenum- 
ber that appears in the governing Helmholtz equation, 



(V 2 + k 2 



u 



(») 



(x) = 0. Another immediate consequence 



is that the intensity correlator is 



»W (x)V n) 



(x + r) 2 \ = A* [1 + 2 4 (kr)] 



(1) 



Berry established the conjecture for an asymptotic 
regime in which wavelengths are much less than system 
size. It is only approximate at finite wavelength, i.e, in 
practice. The conjecture has been shown to be incor- 
rect at finite wavelength, inasmuch as many modes show 
evidence of scarring [jj, existence of which may be un- 
derstood in the light of the "quantum equidistribution 
theorem" put forth by Shnirelman £|. The conjecture 
is manifestly incorrect if attention is restricted to points 
near a boundary where, locally, plane waves are corre- 
lated with their reflections. Nevertheless, numerical and 
experimental evidence shows that it is widely satisfied 
~>\. and references in Q. 

Recent measurements on elastic waves in plates have 
underlined the inadequacy of the conjecture, as stated, 
for systems more complicated than the scalar billiard Q. 
Three dimensional microwave billiards (or merely thick 
quasi 2D billiards) for which the electric field satisfies a 
vector wave equation, and elastic wave systems in gen- 
eral, require a statement about the correlations of the 
vector- valued eigenmodes @. Even in scalar wave sys- 
tems, if they are inhomogeneous, or if interest includes 



points near boundaries, the conjecture needs modifica- 
tion. It is not sufficient to extend the conjecture by 
expressing the modes as uncorrelated superpositions of 
plane waves of all wave types, as their relative ampli- 
tudes remain unspecified. In those recent measurements 
on elastic waves in plates, the intensity correlator did 
not satisfy JQ. This note is intended to provide the ap- 
propriate generalization of Berry's conjecture needed for 
experiments in wave-bearing systems more complex than 
scalar billiards. 

We begin with an identity, written in the form it takes 
for a tensor Green's function appropriate for a vector 
wave equation, 



G (xi,x 2 ,w) = 



E 



u 



M (xi) ®uW (x 2 ) 



(2) 
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as a modal sum over the normalized real modes with 
eigenfrequencies u) n . The imaginary part of this Green's 
dyadic is 

3G ( Xl! X 2 , W ) = ^ Yl U(M) ( X l) ® U( " } ( X =0 6 ^ ~ ^) ■ 



This may be averaged, either over a short range in fre- 
quency, or over an ensemble of systems that differ from 
the system of interest only at positions far from the 
closely spaced points xi and X2. In either case G is 
largely unaffected. The right side becomes the corre- 
sponding modal correlator. It is seen to be in general not 
Jo, but rather 3G. It is not a spatial average, as called 
for by the Berry conjecture, but rather a frequency or 
ensemble average. 

Thus we are led to a generalized Berry conjecture. 
Based upon the exact identity for ensemble or frequency 
averages, we conjecture that it is also true for spatial 
averages at a fixed mode. This reduces to the Berry 
conjecture for the simple case of a scalar wave. The con- 
jecture about the correlator is, as demonstrated above, 
manifestly correct if what one means by the averaging is 
a frequency or ensemble average. The potentially more 
problematic aspects lie in the supposition that this cor- 
relator may be found within spatial averages on a single 
mode of a single sample from the ensemble, or for that 
matter that the statistics are Gaussian. 
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When this generalized Berry conjecture is applied to 
the modes of an infinite isotropic homogeneous 3D elas- 
tic body, the modal correlator is given by the Green's 
function, G°°, satisfying Navier's equations 



(A + /i)V(V-G 00 )+/iV 2 G c 



P uj 2 G c 



-L5 3 (r), 



with I being the identity tensor, A and p being elas- 
tic Lame constants, and p representing material density. 
The spatial Fourier transform of the solution is obtained 
readily: 



pG°° (q,«0 



q® q/ |q|' 



i-q^q/|q| 

|q| 2 c f - ( w + ie f H 2 c 2 t - (uj + ie) 2 



The first term is the longitudinal part with wavespeed 
q, the second is the transverse part with wavespeed c*. 
On taking the imaginary part of the 3D inverse Fourier 
transform, one finds 



3G°° oc J cos (q • r) [c^ 3 6 (|q| - u/d) q ® q/ |q| 2 



c t d <y(|q| -u/a) I-q®q/|q 



^ilqdq. 



It is seen that 3G°° is a superposition of plane waves of 
the two types and of all directions of propagation, with 
relative strengths given by the inverse cubes of the wave 
speeds, i.e. by equipartition (Tll |. 

If the frequency averaging is over a sufficiently broad 
band, then, even in a finite system, the correlator reduces 
to that in the unbounded medium, 3G°°. This is readily 
established, as in 0], by recognizing that short time re- 
sponses are independent of distant parts of a structure, 
and that sufficiently short time responses are equivalent 
to frequency averaging over bands of sufficient width. 
The theorem is readily generalized to the vicinity of a 
boundary, or a scatterer. In these cases it is not difficult 
to show that the diffuse field-field correlator may be also 
constructed by superposing an equipartitioned set of un- 
correlated incident plane or incoming or standing waves 
together with their coherent reflections and scatterings. 
For the special case of elastic waves near a free surface, 
this was explored in (T^ |. 

The Green's function is more complicated in a plate, in 
particular if its thickness is comparable to a wavelength. 
In the work reported by Schaadt et al. Q, an intensity 
correlator was constructed by averages over space and 
over a small number of modes. Due to the good preser- 
vation of up/down symmetry, all modes were either of 
a purely flexural (odd up/down parity) character, or a 
mixture of extensional and shear (even up/down parity). 
Thus their correlator should be the imaginary parts of 
the partial Green's functions. Except for the effects of 
nonzero thickness to wavelength ratio, these modes have 
displacements that are purely out-of-plane or purely in- 
plane respectively. 

Modes which are antisymmetric on reflection about the 
mid-plane, sometimes called flexural, are un-coupled to 



the others; at frequencies below the first cutoff (to — 
nc t /2h, where h is the half- width of the plate) there is a 
single wave number fc/ governing such waves. Their in- 
tensity correlator must therefore be formed from a single 
Bessel function Jo (fc/r). This was in fact observed []]]; 
Schaadt et al. reported a good fit to 1 + 2Jq (fc/r), These 
modes in fact have vector valued fields, so the correlator 
is technically a tensor. The measured correlator is a con- 
traction of that tensor with the (unknown) polarization 
vector p of the detector. Given a field (an eigenmode) 
ip (x) they constructed an 'intensity correlator' 



I (r) = ([p • i, (x)] 2 [p • i, (x + r)] 2 ) / ([ P • i, (x)] 2 )' 



where the overbar indicates additional average over the 
direction of the vector r. We presume that the detector 
polarization p is held fixed during this averaging. Inas- 
much as u is a Gaussian process, the above fourth order 
statistic reduces to the sum of three products of two sec- 
ond order statistics. In terms of the generalized Berry 
conjecture the correlator is rewritten as 



J(r) = l+2(p-3G (x,x + r) • p) 2 / (p • 3G (x,x) • p) 2 . 

(3) 

For the surface of a plate in flexure, even far from the 
edges, such that G « G°°, the Green's function is not as 
simple as might have been hoped, 



3G°° (r) cxxi ® ynv 2 [J (fe/r) - J 2 (fc/r)] /2 



+x 2 ® x 2 ^ 2 [J (fc/r) + J 2 {k f r)\ /2 
+ (xi ® x 3 - x 3 ® xi) vJi (fc/r) 
+x 3 (g) x 3 J (fc/r) . 



Unit vectors Xi and x 2 lie in the mid-plane of the plate, 
with Xi taken in the direction of r, and x 3 is normal to the 
plane and pointing towards the surface at hand. Factor 
v represents a degree of in-plane surface motion associ- 
ated with such waves, and vanishes at long wavelength. 
The correlator I is not only dependent upon separation 
distance r, but also upon the angle between polarization 
of the detector p and the separation vector direction r/r. 

By averaging over the direction of vector r, for a single 
flexural mode, one finds (see Appendix) 



/(r) =1 + 2 



J 2 (fc/r) (1 + ^V/2) 2 + J 2 (fc/r) 

(l + i,V/2) 2 

(4) 

with p 2 = (p 2 +Py) /p 2 - The correlator is plotted in 
Figure for a number of values of vp. 
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Figure 1: Intensity correlator of a single flexural mode for vp 
equal (a) 0, (b) 2, and (c) oo. 

It is equal to 3 at zero separation, r = 0, as de- 
manded by Gaussian mode statistics, and is higher than 
1 + 2Jq (fc/r) for non-zero values of vp, with the most 
pronounced difference observed near the first minimum. 

At realistic values of vp (Schaadt et al. estimate 
p ~ 0.33; we calculate v — —0.68 at the relevant frequen- 
cies), the effect is small, and difficult to resolve within 
the data's precision. The sole anomaly in the data is 
the best-fit value of the relative variance at zero separa- 
tion: 2.93 ± 0.05, an anomaly with only small statistical 
significance. Such a value cannot be explained with the 
current theory; indeed the basic assumption of Gaussian 
statistics demands that this quantity be 3. However, if 
the quantity 2.93 is understood as the ratio between the 
relative variance at zero separation and at the first min- 



imum, then the current theory can explain the anomaly, 
by calling for \vp\ = 1.06. 

The data's precision does not support any more de- 
tailed comparisons. This is also the case with the in- 
plane modes. Modes which have even up/down parity 
consist of an equipartitioned diffuse mixture of longitudi- 
nal waves with wavenumber ki (which have both in-plane 
and, due to Poisson effect, out-of-plane components of 
displacement), and in-plane horizontally polarized shear 
waves with wavenumber k s h- These waves mode con- 
vert to one and other at the plate boundaries. Thus the 
relevant Green's tensor has two wavenumbers, and one 
anticipates structures like those seen in Figs. 4 and 5 
of Schaadt et al. However, as in the flexural case, one 
does not expect to see simple Bessel functions Jo, but 
rather also terms in J\ and J2. The relative amplitudes 
of these several terms are not obvious a priori but could 
be predicted by the present theory. An attempt to fit 
their data to the present theory is probably unwarranted 
at this time. A revisit to structures like theirs, but with 
a well characterized detector of known polarization, may 
be indicated. 

In summary, we have advanced a modification of Berry 
Conjecture, appropriate for the eigenmode statistics of 
wave-bearing systems. It is expected to be relevant, 
not only for elastic waves in homogeneous plates, but 
in general statistical physics of waves in heterogeneous 
and mode-converting systems as well. 
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Appendix A: MULTI-MODE INTENSITY 
CORRELATOR IN A CHAOTIC PLATE 

We start calculation of the full intensity correlator by 
first considering the normal modes of the Rayleigh-Lamb 
spectrum ^3] • The displacement vector of these modes 
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is given by 

u = [U (x 3 ) (k-'V) + x 3 W (x 3 )] f ( Xl ,x 2 ) , 

with / satisfying a scalar 2D Helmholtz equation: 
[V 2 + A; 2 ] f {x\,X2) = 0. The displacement components 
U and W are the solutions of a boundary-value ODE 
in x 3 . With the vertical wavenumbers of longitudi- 



nal and shear waves defined as a 



7< 



k , and 



(3 2 = oj 2 /c 2 — fc 2 , one deduces the dispersion relation for 
the odd and even up/down parity modes 



tan (3h / tan ah = - (fc 2 - 1 ) 2 j Aa(3k 



±1 



where +1 in the exponent corresponds to the odd parity 
modes, and —1 to the even modes. The dispersion re- 
lation gives the wavenumbers of the odd (fc/) and even 
modes (fc;) as multi-branched implicit functions of the 
frequency: k — k n (ui). Expressions for U and W of the 
odd and even modes respectively are 



U ~2k 3 P sin f3h sin ax 3 — (fc 2 — /3 2 ) kj3 sin ah sin j3x 3 , 
W — 2k 2 a(3 sin (3h cos ax 3 + (fc 2 — (3 2 ) fc 2 sinaft,cos/3a;3; 



We construct a partial Green's function (|2J of the 
Rayleigh-Lamb spectrum, and find its imaginary part, 



a n Jo (k n r) 5 a p/2 



+ KJ2 (k n r) (S al3 /2 - rarp/r 2 
3G^=5>„J (fc„r), 

a 

^G^ 3 = - 3G 3q = dn Ji (k n r) r a /r. 



(Al) 



The sum is taken over propagating, i.e. having real 
fc„, modes only. Greek indices span the in-plane 
space: a,/3= {1,2}. By means of the factor, v n = 
U (h) /W (h) |fc=k n (fc,), we can write the modal ampli- 
tudes, a n = b n = c n is 2 7 and d n — c n v n , in terms of 
the amplitude describing out-of-plane displacement of the 
plate surface, 



7r dk k 



4 3uw f +'* 



W 2 (h) 



HP 2 (X3)+W 2 (X 3 )]dx 3 



k—k n (uj) 



U =2k 3 (3 cos /3h cos ax 3 — (fc 2 — /3 2 ) fc/J cos ah cos j3x 3 , 
W — 2k 2 a(3 cos /3h sinaa;3 + (fc 2 — /3 2 ) fc 2 cos ah sin ftx^. 

By specifying a complete set of solutions / in the plane 
(for example, standing plane waves or standing cylindri- 
cal waves), we construct the modes of an infinite plate. 
Alternatively, we may specify a complete set of propa- 
gating waves /, in which case a complex conjugate must 
be inserted on the first factor u in equation QJ. These 
modes are not the natural modes of a finite plate unless 
the boundary conditions at the outer rim are particularly 
special. They may nevertheless be used in a modal ex- 
pansion of the Green's function if attention is confined 
to early enough times (alternatively if a frequency aver- 
aging is done) as discussed above. The average of the 
exact Green's function, G, can then be substituted by 
the Green's function in the infinite plate, G°°. 



The horizontal shear modes have displacements purely 
in the plane of the plate: 

u=^(a;3)(fc- 1 V) x [x 3 f(xi,x 2 )], 

the dispersion relation for the shear wavenumbers (k s h) 

being: fc = fc„ (u) = J (uj/c t ) 2 — (nn/2h) 2 . Imaginary 
part of the corresponding partial Green's function has 
the same form as for the Rayleigh-Lamb modes ijAlll . 
However, the modal amplitudes are now as follows, a n — 
-b n = (1 + S 0n ) /4/ic 2 , and c n = d n = 0. 

The full multi-mode tensor Green's function includes 
the modes of all (namely, odd and even parity Rayleigh- 
Lamb, and horizontal shear) branches required for its 
short-time expansion at a given frequency. The prop- 
agating modes of these branches contribute to the full 
intensity correlator {JJJ, 



I ( r ) = 1-l-2 ( ara< ° 2/2 + Cn) Jo (fc " r)] 2 + E " KJ2 {knT)]2 p4/8 

E„KP 2 /2 + c„)] 2 

I 

The averages over directions of the separation vector r are carried out with the help of the following rules: 

(r a rp/r 2 ) =S a p/2, 
(rarpr^rjr 4 ) = (6a/3&yi + 8 a -ySpi + S al Sp 7 ) /8 
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In the special case that we have, the frequency is such 
that there is only one odd (flexural) mode, and the sum 
is replaced with a single term, yielding correlator 
The factor v is computed for the plate parameters of 



Ref. Q (thickness 3 mm, Poisson ratio 0.33, transverse 
wavespeed, q = 3.1mm//xs), and excitation frequencies 
432, 510, 513, 514kHz, to be v = -0.68. 



